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ABSTRACT 

We present the Morphology-Density and Morphology-Radius relations (T-S and 
T-R, respectively) obtained from the WINGS database of galaxies in nearby clusters. 
Aiming to achieve the best statistics, we exploit the whole sample of galaxies brighter 
than My = —19.5 (5,504 objects), stacking up the 76 clusters of the WINGS survey 
altogether. 

Using this global cluster sample, we find that the T-S relation holds only in the 
inner cluster regions (R < 1/3 x R 200 )) while the T-R relation keeps almost unchanged 
over the whole range of local density. A couple of tests and two sets of numerical 
simulations support the robustness of these results against the effects of the limited 
cluster area coverage of the WINGS imaging. The above mentioned results hold for all 
cluster masses (X-ray luminosity and velocity dispersion) and all galaxy stellar masses 
(M*). The strength of the T-S relation (where present) increases with increasing M*, 
while this effect is not found for the T-R relation. 

Noticeably, the absence/presence of subclustering determines the pres¬ 
ence/absence of the T-S relation outside the inner cluster regions, leading us to the 
general conclusion that the link between morphology and local density is preserved just 
in dynamically evolved regions. We hypothesize that some mechanism of morphological 
broadening/redistribution operates in the intermediate/outer regions of substructured 
(“non relaxed”) clusters, producing a strong weakening of the T-S relation. 

Key words: galaxies: clusters - galaxies:general - galaxies: morphology 


1 INTRODUCTION 

The detection of significant differences in the galaxy 
populations between clust ers and field dates back to 
iHubble and Humasonl rtl93ll ). who noted that “the predom- 

* E-mail:giovanni.fasano@oapd.inaf.it 


inance of early types is a conspicuous feature of clusters in 
general”. Twenty years later, the prevalence in clusters of 
stellar populat ion II galaxies (early type s, mainly SOs) was 
interpreted by ISpitzer and Baade! 11195 lh as due to galaxy 
collisions in the dense cluster environment. Attempts to clas¬ 
sify clusters according to the galaxy morphological composi- 
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tion were subsequently made bv lMorga ^1 (Il96lh and lOemleil 
(Il974h . 

Yet, the first quantitative assessment of the relation be- 
tween galaxy m orphology and environment was provided by 
iDresslerl (Il980l . D80 hereafter), based on photographic plates 
of 55 rich galaxy clusters, obtained using the Las Campanas 
2.5m and the Kitt Peak 4m telescopes. Indeed, it was after 
this classical paper that the relation between the local pro¬ 
jected density of galaxies in the sky (computed using the 10 
nearest neighbours; Eio hereafter) and the fraction of galax¬ 
ies of different morphological types (the morphology-density 
relation; T-E hereafter) began to be considered among the 
most important empirical evidences of extragalactic astron¬ 
omy. In particular, the D80 finding that the T-E relation 
holds in both structurally regular (dynamically “relaxed”?) 
and irregular (clumpy, “non relaxed”) nearby clusters and its 
rapid extention to the group and gene r al field environments 
(Bh avsa 3 Il98ll : [de Souza et alJ 1 19821 ; [Postman and Gellei) 
1 19841) . led the astronomical community to consider this rela¬ 
tion a sort of universal rule in the extragalactic astronomy. 
Indeed, the fact that early-type galaxies (Es and SOs; ETGs 
hereafter) preferentially dominate dense environments, while 
late-type galaxies (Spirals and Irregulars; LTGs hereafter) 
are more common in low-density regions has progressively 
become an out-and-out paradigm with which to reckon when 
dealing with problems related to the evolution of galaxies, 
from both the phenomenological and theoretical side. For 
instance, lEvrard et al] (fl990r) explored a model in which 
galaxy morphology is intrinsically related to the peak height 
of the initial fluctuations in a CDM universe, finding good 
agreement with D80 for the T-E relation of Es. 

Besides the existence of the T-E relation in galaxy 
clusters, in D80 it is also claimed that this relation is 
much stronger than the relation between morphology and 
cluster-centric distance (morphology-radius relation; T-R 
hereafter) . Ho wever, a dif feren t viewpoint was proposed by 
IWhitmore and Gilmorel (Il99ll h who re-examined the D80 
data, finding that the strengths of the T-R and T-E rela¬ 
tions are comparable. They also suggested that the physical 
mechanisms driving the gradients f or Es are dif ferent from 
t hose operating for spirals and SOs. IWhitmore and Gilmorel 
dl993h examined again the D80 sample concluding that the 
T-R, normalized to a characteristic radius (R c ), is the fun¬ 
damental relation driving the morphological fractions of all 
types, from spiral galaxies (slightly decreasing toward the 
center, then going close to zero near the center), to SOs (mod¬ 
erately increasing toward the center for R/R c <0.2, then 
sharply falling down) and Es (from 10% in the outer regions 
to ~16% at R/R c ~1, then rapidly in creasing up to 60-70% 
in the inner regions). Moreover. IWhil.morel (119951) analysed 
compact and loose groups and the field, finding that the T- 
S is not present in these environments, just strongly rising 
near the center of clusters. 

The Whitmore & Gilmore’s viewpoint about the preva¬ 
lence of the T-R over the T-E relation in determining the 
morphological fractions in clusters was almost bailed out 
after the T-E relation was proved by the MORPHS group 
dDressler et al.lll997l . D97 hereafter) to hol d also for interm e- 
diate redshift clusters (z ~0.5; see also lFasano et al . 200(3 . in 
the redshift range 0.1-0.25). In this redshift range, however, 
the relation is found very weak or even absent for irregu¬ 
lar clusters. Moreover, we mention that iTreu et ~ahl (|2003l ) 


claimed the cluster-centric distance to be the driving pa¬ 
rameter of the ETG fraction in the inner part of the cluster 
Cl 0024+16 at 2 ~0.4. 

Although structure and morphology are intrinsically 
different properties, thus making dangerous to compare vi¬ 
sual T-E (morpholog y) at low z with s tructural T-E (Ser- 
sic index) at high 2 ( van der Well 120081 ) . in the subsequent 
decade several authors studied the T-E relation in different 
ranges of redshift and stellar mass. Most of them agreed 
that: ( i ) the T -E is already in plac e, although in a less 
strong fashion dPostman et al.l 120051 ). at 2 ~0. 5, both in 
groups and (possibly fr o m 2 ~1) in clust ers dTreu_gLal 


2003; ISmith et alJ 120051 : iDesai et al.l 120071 : IWilman et al 


20091) ; (ii) the evolution o f morphological fract i ons strongly 


depends on galaxy mass dNuiiten et al.l 120051 : lOesch et al.l 
l20ld : IVulcani et al.l 1201 ll ) ; (iii) in dense environments the 
ETGs dominate a t all redshift s and their fraction in¬ 
creases over time dNuiiten et al.l 120051 : ISmith et al.l 120051 : 
ICapak et al.ll2007 ), mostly because of the grow ing fraction of 
SOs ( Postman et al.ll2005l : IVulcani et al.ll201ll) : (iv) the evo¬ 
lution of the T-S relation turns out to be less st rong for mass 
limited than for magnitude l i mited samples jHolden et al.l 
120071 ; Ivan der Wei et al1l2007l : iTasca et al .II2OO9I + 

Besides the T-E dependence on redshift, some pa¬ 
pers have also analysed the morphological fractions as a 
function of the g lobal cluster properties. In particular, 
IPesai et alJ (12007 ) found a mild dependence of the mor¬ 
phological fractions on the cluster velocity dispersion, while 
the growt h of the Spir al/SO fraction with redshift was 
found by lPoggianti et ahl |2009l ) to be stronger for low mass 
than for high mass clusters. A specific effect of the clus- 
ter environment on the Spiral/SO fraction is claimed by 
ICalvi et al.l (120121 ). based on the sharp enhancement/dearth 
of SOs/late-types found in clusters at decreasing redshift, 
compared to other environments (single galaxies, binary sys¬ 
tems, poor/rich groups). 

A totally different and innovative approach to the T-E 
has been proposed with the kinem atic morphology-density 
relation bv ICappellari et al.l (1201 ll . Paper VII of the Atlas3P 
collaboration). Rather than rely on the visual morphology, 
they classify galaxies according to their stellar kinematics, 
dividing ETGs into slow- and fast-rotators (instead of Es 
and SOs). In this way and adopting a very local density es¬ 
timator, computed just using the 3 nearest neighbours (the 
Atlas3P galaxies mainly inhabit low-density environments), 
they found a T-E cleaner than using classic morphology. 
They interpret this improvement as due to the fact that 
many fast-rotators, visually classified Es are actually face- 
on, misclassified SOs. 

WINGS is a multiwavelength photometric and spec- 
troscopic survey o f 77 galaxy clusters at 0.04<z<0.07 
dFasano et al. I [20061 F06 hereafter). Clusters were selected 
in the X-ray from th e ROSAT Brightes t Clu ster sam¬ 
ple and its extension (lEbeling et all Il998. 200 (3) and the 
X-ray Brightest Abell-type Cluster sample ( Ebeling et al.l 


fBV) for all 77 fields (Varela et al.l 

200(3, WINGS-OPT), as 

well as infrared (JK) photometry 

Valentinuzzi et al. 

2009, 

WINGS-NIR), multi-fiber snectroscoov (Cavaetal. 

2009, 

WINGS-SPE) and U-band photometry (Omizzolo et al. 


face photometry and morphology of WINGS galaxies have 
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been obtained thr ough the purposely d evised, automatic 
tools GASPHOT (ID Onofrio et alJlioTil and MORPHOT 
llFasano et al.l 12012 . F12 hereafter), respectiv ely. A full de¬ 
script ion of the WINGS database is given in lMoretti et all 
ll2014ll . 

Exploiting the large and coherent database provided by 
the WINGS project, we revisit here the T -E and T-R rela¬ 
tions in galaxy clusters of the local Universe. In particular, 
thanks to the large galaxy sample, we are in the position 
of testing how these relations behave in different ranges of 
cluster-centric distance (for T-E) and local density (for T- 
R), as well as for galaxies in different stellar mass ranges 
and for different global cluster properties. In Section [2] we 
describe the cluster and galaxy samples. In Section [3] we 
give detailed information about the data used to perform 
the analysis. The procedure used to account for and remove 
field galaxies is described in Section [4] Section [5] and Ap¬ 
pendix are devoted to characterize the issue related to 
the limited cluster area coverage of the WINGS imaging 
and to investigate the effect of such an issue on the results 
we obtain for the T-E and T-R relations. These results are 
presented in Sections[6] In particular, in Sections I6.3l and r6.4l 
we investigate how T-E and T-R depend on the global clus¬ 
ter properties and on the galaxy stellar mass, respectively. 
In Section [7] we summarize our results and try to outline 
possible scenarios able to interpret them. 

Throughout the paper we adopt a flat geome¬ 
try of the Universe with the following cosmology: 
Ho=70 kms _1 Mpc _1 , fiM=0.3 and flA=0.7. 

2 THE SAMPLE 
2.1 The cluster sample 


we assume that a failure of the GASPHOT fitting occurred, 
thus preventing us from relying on the output GASPHOT 
parameters. For the resulting sample of 39,3 04 o b jects, 
putti ng together the data from WINGS-SPE (tCava et al.l 
1200914 and from the literature, we gathered spectroscopic red- 
shifts of 9,361 galaxies, of which 6,297 are cluster members. 
The absolute magnitudes of all galaxies in our sample, but 
the spectroscopic non-members, are computed usin g the av¬ 
erage redshifts of clusters given in I Cava et al.l (120091 4 and the 
cosmology gi ven in Sect ion [Tl They are K-corrected using 
the tables in IPoggiantTI (ll997t) and according to the mor¬ 
phological types provided by MORPHOT (see Section mi) . 

When comparing our results with those reported in the 
literature, we will mainly advert to D80 and D97, which 
still remain the most important and cited references on this 
subject. After accounting for the different waveband and 
cosmology, the absolute V-band magnitude limits (My m ) 
given in these papers for the galaxy samples are —19.67 
and —19.17, respectively. We adopt M y m = —19.5, which 
is also close to the value u sed for the ENACS sample by 
iThomas and Katgertl d2006l . ;M y rn = —19.57). In this way, 
we are left with the final ’reference sample’ of 5,504 galaxies, 
of which 82% have measured redshift and 67% are cluster 
members. 

Throughout the paper, alternative sample partitions 
and selection criteria will be tested, so that the actual sam¬ 
ple size might vary from time to time, also depending on 
the expected field counts relative to the particular selec¬ 
tion. For this reason, we usually report each time in the 
following tables the proper sample sizes before the field re¬ 
moval. Finally, the same above mentioned limit of luminos¬ 
ity (My™ = —19.5) is adopted for neighbour galaxies in the 
computation of the local densities ISection 13.21) . 


The galaxies used for the present analysis are located in the 
fields of 75 clusters belonging to the WINGS-OPT survey 
(see Table 5 in F06). The cluster A3164 has been excluded, 
due to the poor PSF quality which prevented us from obtain¬ 
ing a reliable morphological classification. Moreover, when 
performing the analysis of the T-R relation, we excluded the 
galaxies in the field of clusters A311, A2665 and ZwC11261. 
In fact, for these clusters we do not have a redshift cover¬ 
age sufficient to allow a reliable estimation of the velocity 
dispersion, which is needed to compute the characteristic 
cluster radius R 200 (as defined in Sec. 1531 and normalize 
the distance of galaxies from the cluster center. 

2.2 The galaxy sample 

Our starting sample of galaxies roughly concides with that 
illustrated in F12. More precisely, the MORPHOT cata¬ 
log of 39,923 WINGS galaxies has been cross-matched with 
the catalog of structural parameters obtain ed running GAS- 
PHOT onto the V-band WINGS imaging llD’Onofrio et al.l 
l2014l l. This catalog provides apparent total magnitudes and 
structural parameters of galaxies. The apparent magnitudes 
ha ve been correc t ed for galactic extinction using the maps 
in ISchlegel et all lll998h . The galaxies for which the Sersic 
index derived by GASPHOT turned out to coincide with 
the boundary values of the allowed search interval (0.5-8) 
have been excluded from the database. In fact, in these cases 


3 THE DATA 

All quantities we describe in the followings sub-secions and 
we use in the following analysis can be easily retrieved, 
through the Virtual O bservatory tools, fol lowing the rec¬ 
ommendations given in lMoretti et al.l (120141 4. 

3.1 Morphology 

The morphological types of galaxies are taken from the 
database of WINGS galaxies at: 

http://cdsarc.u-strasbg.fr/viz-bin/VizieR?-meta.foot 
&-source=J/MNRAS/420/926. This database has been ob¬ 
tained running MORPHOT (F12) on the V-band WINGS- 
OPT imaging of galaxies with isophotal area larger than 
200 pixels at the threshold of 2.5<7bk g (where tJbkg is the 
standard deviation of the background). It is worth noticing 
that none of the WINGS-OPT galaxies with My ^ —19.5 
has been removed from the final sample because of the 
200 pixels cutoff adopted by MORPHOT. This means that 
our absolute magnitude limited sample is safe from surface 
brightness incompleteness (very compact galaxies). 

MORPHOT is an automatic tool purposely devised 
to obtain morphological type estimates of galaxies in the 
WINGS survey. It combines a large set of diagnostics of 
morphology, easily computable from the digital cutouts of 
galaxies. MORPHOT produces two independent estimates 
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of the morphological type based on: (i) a semi-analytical 
Maximum Likelihood technique; (ii) a Neural Network ma¬ 
chine. The final estimator has proven to be almost as ef¬ 
fective as the visual classification. In particular, it has been 
shown to be able to distinguish between ellipticals and SO 
galaxies with unprecedented accuracy (see Fig. 11 in F12). 

The MORPHOT estimator (MType in the database) 
is almost coincident with the Revised Hubble Type defined 
bv lde Vaucouleurs et al.l dl99ll . see Table 1 in F12). For the 
purposes in the present paper it is convenient to split the 
galaxy sample in three broad morphological types: (i) ellip¬ 
ticals (E); (ii) lenticulars (SO); (iii) spirals (Sp). With this 
convention, the (few) irregular galaxies in our sample are 
included in the Sp class and the cD galaxies are included in 
the E class. 

Using a sample of 176 common galaxies belonging to 
18 common clusters, we have compared the broad morpho¬ 
logical type (Es/SOs/Spirals) provided by MORPHOT with 
the corresponding classifications given in D80. The agree¬ 
ment turned out to be quite good for spirals (discrepancy of 
2%), while for Es and SOs we found a discrepancy of ~20%, 
mainly in the sense of an excess of elliptical galaxies in the 
MORPHOT classifications with respect to the D80 ones. 
The most obvious explanation of such discrepancy is that it 
originates from the different material used to inspect clus¬ 
ters: photographic plates in the original paper of D80 and 
CCD imaging in the WINGS survey. For this reason, one 
of us (AD, indeed!) has visually (and blindly) re-classified 
25 discordant objects using the WINGS observing material 
(CCD). With the new AD classifications the discrepancy 
turned out to be reduced down to a physiological amount 
(~6%)- 


3.2 Local Density 

The projected local galaxy density relative to a given galaxy 
is commonly defined as the number of neighbours (Nn)of 
the galaxy per square megaparsec. According to the pre¬ 
scriptions given in D80, we compute the projected local 
densities using the 10 nearest neighbours of the galaxy 
with Mv ^ —19.5. At variance with D80, we use a cir¬ 
cle (rather than a rectangle) to compute the area including 
the 10 nearest neighbours. That is: Eio = 10/Aio, where 
Aio = nR 2 0 (Mpc) and Rio(Mpc) is the radius (in Mpc) of 
the smallest circle centered on the galaxy and including the 
10 nearest neighbours. Numerical simulations have shown 
that using the circular area produces a negligible (~0.03) 
downward shift of LogEio with respect to the values com¬ 
puted using the rectangular area. 

In order to perform the correction for field contamina¬ 
tion, rather than our incomplete membership information, 
we prefer to use a statisti cal co rrecti on based on the field 
counts provided bv lBerta et al.l (l2006l f using the deep ESO- 
Spitzer wide-area Imaging Survey (ESIS). In particular, we 
compute for each cluster the number of field galaxies per 
Mpc 2 (AV) up to the apparent magnitude corresponding to 
My = —19.5 at the cluster redshift. Moreover, in counting 
the neighbours within a circle centered on a given galaxy, we 
statistically account for the objects possibly lost because of 
the proximity of the galaxy to the edges of the fied of view 
(including the gaps between the chips). We do this dividing 



Figure 1. Difference between the Eio values obtained with and 
without SDSS information as a function of the coverage fraction. 

the counts by the ratio f c (^1) between the area actually 
covered by the detector and the circular area. 

Starting from the nearest neighbour (i=l) and progres¬ 
sively including the next nearest one (i=i+l), we compute 
each time the circular area A; (in Mpc 2 ) relative to the 
neighbour distance and the corresponding number of neigh¬ 
bours: Ni : n = i/f c — Nf x Ai. When Ni }Tl overcomes 10, the 
final area Aio is computed interpolating between Ai and 
Ai- 1 - 

For 24 clusters in common with SDSS-DR7, we com¬ 
plement our positional and photometric database with 
SDSS in format ion. In particular, we use the relations in 
iFukugita et ali (l2007fl to get the V-band magnitudes of 
SDSS galaxies from the g' and r' ones. For this galaxy sam¬ 
ple we are allowed to obtain complete coverages within Aio 
even for galaxies close to the frame borders, thus being able 
to compute the Eio values (E sdss) without correcting for 
the coverage fraction (/ c =1). Figure [l] illustrates how the 
difference between the Eio values obtained with and with¬ 
out SDSS information depend on the coverage fraction. It 
is evident from Fig. |T]that for small values of f c the statis¬ 
tical correction tends to produce an overestimation of Eio. 
This is likely due to the negative outward density gradient 
in clusters, which makes the real number of galaxies lost be¬ 
cause of the proximity to the frame edges, lower than the 
number obtained correcting with f c , i.e. assuming a uniform 
density inside the circles. For this reason the local densities 
of galaxies for which no SDSS information is available, have 
been statistically corrected using the best-fit relation illus¬ 
trated in Fig. [T] LogEio rr '=LogE 10 +0.178 (f c - 1). 

3.3 Cluster-Centric Distance 

We derive two different cluster-centric distances ( R ), de¬ 
pending on the choice of the cluster center. The first one 
(Rb) assumes the brightest cluster galaxy (BCG) as clus¬ 
ter center, while the second one assumes the center of the 
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Figure 2. Map of the cluster Abell 1668, obtained using the photometric data from the WINGS-OPT catalog and the morphologies 
provided by MORPHOT. The different broad morphological classifications are reported with different colors: red for Es, green fro SOs, 
cyan for early spirals, blue for late spirals and irregular galaxies. Small crosses refer to galaxies fainter than the limit we adopt in our 
analysis (My > —19.5, adopting the cluster redshift distance). Galaxies brighter than this limit are represented by ellipses, whose shapes 
obey the structural parameters provided by SExtractor (isophotal area, axial ratio and position angle). Spectroscopically confirmed cluster 
members and galaxies without spectroscopic information are indicated by fu ll and empty ellipses, respectively, while black symbols (both 
ellipses and crosses) identify galaxies which, according to ICava et al.l i2009h , are considered spectroscopic non members. Finally, dashed 
circles correspond to different distances from the BCG: 0.33, 0.5 and 1 (in R 200 units). 


cluster to coincide with the maximum intensity of the X- 
ray emission. In any case, the cluster-centric distances have 
to be normalized to some characteristic cluster size. One 
class of size estimates is based on top-hat filtered spher¬ 
ical overdensities. In this model, clusters are expected to 
be virialized within regions where the enclosed me an mass 
density exceeds the critical density by a factor 200 (I Peebles! 
ai eacoc 3 119991 , p. 25 and 15, respectively). We assume 
the radius at which this overdensity is reached, R 200 , as the 
characteristic radius of the cluster. The R 200 values for our 
clusters are compute d from th e WINGS cluster velocity dis- 
persions (a) give n in ICava et ahi (I2009I) . using the formula 
jFinn et alJl200j ): 


R 200 = 1-73 x 


1000 


kmS ^’ y/^A + lWl +W 


h - 1 (Mpc) 


Therefore, in each cluster, the cluster-centric distances 
of galaxies (R), to be used in the following analyses, have 
been obtained dividing the distance in Mpc by the proper 
value of R 200 ■ 


I 11 Figure [2] we illustrate, as an example, the map of 
the cluster Abell 1668, obtained using the photometric data 
from the WINGS-OPT catalog and the morphologies pro¬ 
vided by MORPHOT (see the figure caption for details 
about symbols). Similar maps of all WINGS clusters (76) 
are available in the online material (see the electronic ver¬ 
sion of the paper). 


4 REMOVING THE FIELD GALAXIES 

In order to obtain reliable T-E and T-R relations in nearby 
clusters, it is important to properly remove the contribution 
of field galaxies (mostly background) in each broad mor¬ 
phological class (E/SO/Sp). The percentages of field galax¬ 
ies in each class have been estimated u sing the Pa dova- 
Mille nnium Galaxy and Group Catalogue llCalvi et al.ll20lil . 
12011 PM2GC), consisting of a spectroscopically complete 
sample of galaxies at 0.03^ z ^0.11 brighter than Mb = 
— 18.7. Thi s sample is sourced from the Millenn ium Galaxy 
Catalogue jLiske et al JI 200 I iDriver et al.ll2005l ). a B-band 
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Figure 3. Morphological fractions of PM2GC galaxies as a func¬ 
tion of the apparent V-band magnitude. Error bars report the 



Figure 4. Counts of field galaxies from [Berta et aD (;200(> ) and 
from PM2GC. The vertical dotted line corresponds to the ap¬ 
parent magnitude at which we have computed the morphological 
fractions (V=17.77). 


contiguous equatorial survey, complete down to B=20 and 
representative of the general field population in the local 
Universe. The advantage of using PM2GC data is that both 
the instrumental set up of the imaging (WFC@INT) and the 
tool used to estimate the morphological types (MORPHOT) 
are the same for both PM2GC and WINGS galaxies. This 
guarantees a full morphological consistency between the two 
samples. 

Figure [3] reports the morphological fractions of the 
PM2GC galaxy sample in different bins of apparent mag¬ 


nitude. Even if the expected increase of the E fraction and 
decrease of the Sp fractions are clearly detectable, Fig. [3] 
shows that the percentages E/SO/Sp depends weakly on the 
apparent magnitude. In our analysis we assume E/SO/Sp 
= 19.3(±3.1)%/21.7(±2.7)%/59(±5.3)%. These values cor¬ 
respond to the morphological fractions of PM2GC galaxies 
with V^17.77, this limit including 95% of the WINGS galax¬ 
ies in our reference sample (see below). 

In Figure [4] the global field counts in the V-band from 
the PM 2GC galaxy sam ple are compared with the counts 
given bv lBerta et alj (120061 1 in the same waveband and in the 
magnitude range of interest for us. The slight systematic ex¬ 
cess of the P M2GC counts for bright galaxies with respect to 
the counts bv lBerta et al.l ()2006l l is likely due to the presence 
in the PM2GC field of many galaxy groups, which, because 
of the criteria adopted in the choice of th e fie ld, are proba¬ 
bly lacking in the ESIS database used bv lBerta et al.l (|2006T) . 
On the other hand, the sharp decline of the PM2GC counts 
towards faint galaxies clearly reveals the incompleteness of 
the PM2GC galaxy sample for which a morphological clas¬ 
sification has been possible through MORPHOT (isophotal 
area larger than 200 pixels at the threshold of 2.5abk g )- Nev¬ 
ertheless, it is noticeable that the two counts curves assume 
roughly the same value (1.51, that is ~32 galaxies per square 
degree) at V^17.77, which is the value chosen above to com¬ 
pute our morphological fractions. All things considered, the 
statistical counts of field galaxies suggest us to remove in our 
analyses 32 galaxies per square degree, in particular: 6 Es 
(19.3%), 7 SOs (21.7%) and 19 Spirals (59%). 

According to Section 12.21 in our galaxy sample the 
global spectroscopic coverage and cluster membership, up 
to the adopted limit of absolute V-band magnitude (—19.5), 
are 82% and 67%, respectively. These rather high percent¬ 
ages allowed us to deal with the problem of field galaxies 
using a strategy complementing the statistical field counts 
with the membership information. In particular, we first re¬ 
moved from the sample those galaxies which, according to 
the redsh ift in format ion and the membership criteria de¬ 
scribed in lCava et al.l (l2009l l. are not cluster members. Then, 
if the number of galaxies removed from each broad morpho¬ 
logical class (Es/SOs/Sp) turned out to be lower than the 
expected field counts relative to that class, we further sub¬ 
tracted from the residual t he quantity needed to reach the 
field counts expected from IBerta et al.l |200d l. In this way, 
it might happen that the actual number of removed galaxies 
exceeds the number expected from the field counts. In any 
case, we have verified that the differences between the T- 
E relations obtained using our mixed strategy and a ‘pure’ 
statistical approach are practically negligible. 


5 THE CLUSTER AREA COVERAGE ISSUE 

Having at our disposal a large galaxy sample in a well de¬ 
fined, complete cluster sample, our strength lies in the op¬ 
portunity to take advantage of a good statistics. For this 
reason, in the present analysis we mostly use galaxies in the 
whole cluster sample. At times, however, in order to per¬ 
form some particular analysis, we use galaxies belonging to 
a sub-sample of clusters (selected, for instance, on the basis 
of X-ray emission, velocity dispersion or subclustering), or 


© 2012 RAS, MNRAS 000, []]-?? 
































1 


Morphological Fractions in WINGS 7 



0.35 0.5 0.75 1 1.35 1.5 

R b (distance from BCG in R 200 units) 


Figure 5. Upper Panel: The fraction of WINGS clusters with 
circular coverage fraction greater than some given value (0.50, 
0.75, 0.90 and 0.95) as a function of the radius (in units of 
11/It 200 ) of the circular area centered on the BCG. Bottom 
Panel: The solid line reports the fraction of WINGS clusters 
whose equivalent radius (see text) exceeds the value of Rg re¬ 
ported in the abscissa. The dashed line (red in the electronic 
version) illustrates the fraction of WINGS clusters for which 
the largest clustercentric galaxy distance ( Rg ) exceeds the value 
given in the abscissa. 


galaxies obeying some constraint (for instance on the stellar 
mass). 

On the other hand, our weakness lies in the limited and 
not always regular cluster area covered by the WINGS im¬ 
ages. The upper panel of Figure [5] shows that, for Rs=0.5, 
0.75 and 1.0, the WINGS images cover more than 90% of 
the circular (BCG centered) area in just ~66%, 28% and 
2% of the clusters, respectively, while a coverage fraction 
(CF hereafter) of 0.75 at the same radii is reached in ~85%, 
50% and 11% of the clusters. Similar CFs are found if we 
assume the cluster centers to coincide with the maximum 
intensity of the X-ray emission. An alternative way to illus¬ 
trate the coverage problem is provided by the bottom panel 
of Fig. [5] In this panel the solid line reports the fra ction of 
WINGS clusters whose equivalent radius R e , = tUoo/tt 
(A 200 being the area covered by the WINGS images in units 
of R 200 ) exceedes the value of Rg reported in the abscissa, 
while the dashed line (red in the electronic version) reports 
the fraction of WINGS clusters for which the galaxy with 
the largest clustercentric distance has a Rg exceeding the 
value given in the abscissa. 

Even if one of the goals of the present paper is the analy¬ 
sis of the morphological fractions as a function of the cluster¬ 
centric distrance ( T-R ), from previous studies we know that 
the spiral fraction inside a given (BGC or X-ray centered) 
cluster area is an increasing function of the radius. Thus, 
the limited cluster area coverage of the WINGS images in 
comparison with other analyses (e.g.: D80: lGoto et al.ll2003l : 


iThomas and Katgertl 120061 ). is expected to result in a re¬ 
duced global fraction of spirals. 

Besides this obvious effect, additional biases due to the 
limited cluster area coverage of our data could affect the 
T-R and T-X relations we present here. 

Concerning the first relation (T-R), there is no reason to 
believe that any galaxy included in our sample, with cluster¬ 
centric distance (R) not completely covered by our imaging, 
should be drawn from a morphological distribution differ¬ 
ent from the average distribution at that R. Therefore, the 
T-R should be unaffected by both the limited cluster area 
coverage of our images, and their irregular shape (for INT 
imaging, in particular). For the T-R, the above issues should 
just translate into a limited extension of the relation and a 
relatively large poissonian uncertainty of its outer part. 

Trying to quantify the possible biases affecting the T- X 
relation because of the limited cluster area coverage of our 
imaging, in Appendix[A]we present three different tests. The 
first two tests investigate the effects on T-X of two different 
issues: the irregular shape of WINGS images and their lim¬ 
ited coverage of the cluster area. The third test, performed 
through quite simple numerical simulations, is again aimed 
at investigating the effect of the limited cluster coverage on 
the T-X in different ranges of clustercentric distance. More¬ 
over, since the simulations require some hypothesis about 
the driving parameter (X or R) of the morphological frac¬ 
tions, they could also provide useful indication with regards 
to such parameter (see Section fA3l in Appendix EJ). 


6 RESULTS: THE T-X AND T-R RELATIONS 
IN WINGS CLUSTERS 


In Fig. [o] we have shown that, due to the limited field of 
view of the cameras used in the WINGS optical survey, 
the cluster area coverage of our cluster images turns out 
to be relatively small. As expected, this makes the global 
morphological content of our galaxy sample quite different 
from that reported in previous analyses which made use of 
much larger field of view imaging. In particular, we find 
E/S 0/Sp~33%/44%/23 % (see also IPoggianti et al.l I2009I . 
and lVulcani et aJJ 12011. with a slightly different sample se¬ 
lection), to be compared with the morp hological fractions 
given by D80 (~18%/41%/41%) and IThomas and Katgertl 
112006I . ~19%/47%/34%). It is worth stressing that, sup¬ 
ported by the comparison between the MORPHOT and D80 
morphological classifications (see Section 13.111 . we entirely 
attribute to the limited cluster area coverage of the WINGS 
imaging the above discrepancy of global morphological frac¬ 
tions. 

Let’s now indicate with FfE/Ffso/Ffs P =0.193/0.'217/0.59 
the morphological fractions we have previously found for 
field galaxies, and with Ntot and Nf the total number of 
galaxies and of field galaxies relative to the particular selec¬ 
tion (of clusters and/or galaxies and/or range of I? or X 10 ) 
we are testing in the analysis of the T-X and T-R relations. 
For each bin ( i ) of LogXio or Rg and for each broad morpho¬ 
logical type (T=E/S0/Sp), we compute the morphological 
fractions as: 


Ft = 


(1) 


In these expressions, the symbols N and F indicate 
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Log (S,„) 

Figure 6. T -S (upper panel) and histogram of local density 
(lower panel) for the whole ’reference sample’ of 5504 WINGS 
galaxies with My < —19.5. Full dots, triangles and squares (red, 
green and blue in the electronic version of the paper) refer to Es, 
SOs and Spirals, respectively. The error bars correspond to poisso- 
nian uncertainties. Details about the coefficients reported in the 
upper panel (KS, ACC and ASlope) are given in the text. 

galaxy numbers and fractions, respectively. In particular, 
N l and are the total number of galaxies and the num¬ 
ber of galaxies with broad morphological type T in the i th 
bin, while we indicate with Nf the product N' x Ff and 
with Ff the fraction Nf/Ntot ■ In the above formulation of 
we assume that the number of field galaxies per square 
degree and the global morphological field ratios F/t do not 
depend on the local density. In the following figures we just 
plot the bins where the number of cluster galaxies (N 1 — Nf) 
is greater than 10. 


6.1 T- E relation 

Figure [6] illustrates the T- E relation for the whole reference 
sample of 5,504 WINGS galaxies (see also Table [T|), adopt¬ 
ing a bin size of 0.2 in LogEio. In the upper panel of the 
figure and in the following T- E (and T-R) plots, we report a 
few coefficients, through which we try to quantify the differ¬ 
ences among the Fraction-LogEio (Fraction- Rb for the T- 
R) relations relative to Es, SOs and Spirals. In particular, for 
each pair of broad morphological types, we report the two- 
samples Kolmogorov-Smirnov probability (2S-KS hereafter) 
that the two distributions of LogEio are drawn from the 
same parent population. Moreover, we compare the weighted 
T-E of Es and Spirals, reporting the coefficients ACC and 
ASlope. They give respectively the halved difference of Cor¬ 
relation Coefficients and the Slope difference between the 
T-E of Es and Spirals assuming a linear fit. These values, 
together with their expected r.m.s. uncertainties, are also 
reported in Table [T| 

Figure [S] clearly shows that for the WINGS cluster 
galaxies we recover the classical T-E relation: at increasing 
local density, the fractions of Es and Spirals increase and 
decrease, respectively. The correlations are very strong in 
both cases, as indicated by the coefficients ACC and ASlope, 
while the KS(E-Sp) probability indicates that Es and Spirals 



Figure 8. T-R (upper panel) and histogram of Rg (lower panel) 
for the whole sample of 5187 WINGS galaxies with My <1 —19.5 
for which the value of Rg in units of R 200 is available. The mean¬ 
ing of the symbols is as in Fig. [TT| 

are quite distinct from each other. It is also worth noticing 
that the fraction of SOs seems not to depend at all on the 
local density, as found also by D97. In addition, the KS 
analysis suggests that the SOs constitute a population well 
distinct from both Es and Spirals. The first row of Table [T] 
(columns 2, 3 and 4) reports the values of the coefficients, 
together with their r.m.s uncertainties. 

Let’s now check out whether and how the T-E relation 
depends on the position of the galaxy inside the cluster, in 
particular on the clustercentric distance relative to the BCG 
(Rg). Figure [7] illustrates the T-E in two different ranges of 
Rb. The left panel of the figure refers to the central region 
of the clusters (Rg <0.5), while the right panel refers to the 
Rg interval (0.5-1). It stands out that, while in the central 
part of the cluster the T-E is quite strong, when moving 
outside it seems to become very weak or even absent. This 
is formally confirmed by the ACC and ASlope values, as 
well as by the 2S-KS probabilities obtained in the two Rb 
intervals. These values are reported in Tabic |T] (rows 2 and 
3; columns 2, 3 and 4). 

Tablc[l]also reports, in rows 4, 5 and 6 (columns 2, 3 and 
4), the results we obtain splitting the Rg interval (0-1) into 
three parts, while in the row 7 we report the data relative 
to the Rg interval (0.33-1), which we refer to in Figure [T^J 
Table Q] clearly shows that the region where the T-E actually 
operates is the very inner cluster region (Rg <l/3). The 
high significance found for the T-E for the whole galaxy 
sample (see Fig. © is just due to the fact that the inner 
regions are the most densely populated in the clusters. 

6.2 T-R relation 

In the analysis of the T-R relation we assume again the 
cluster center to coincide with the position of the BCG (R = 
Rb), while in this case we adopt a bin size of 0.15 for Rg. 
We note, however, that both the previous and the following 
results remain practically unchanged if we adopt a different 
bin sizes and/or assume the cluster center to coincide with 
the maximum intensity of the X-ray emission. 
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Log (E 10 ) Log (E, 0 ) 


Figure 7. T-E for WINGS galaxies with Rg <0.5 (in units of R 200 , left panel) and 0.5- Rg <1 (right panel). The meaning of the 
symbols is as in Fig. [6] 


Table 1. Coefficients ACC and ASlope for T-E and T-R in different intervals of Rg and LogEio, respectively. The table also report the 
2S-KS probabilities that the LogEio or Rg distributions of elliptical and spiral galaxies are drawn from the same parent population for 
the T-Y, and T-R relations, respectively. The number of galaxies used in each case is reported below the Rg or LogEio intervals, while 
the expected r.m.s. uncertainties of the coefficients ACC and ASlope are reported below the relative values. 




T-Y 




T-R 


Rg 

ACC 

ASlope 

KS(E-Sp) 

LogEio 

ACC 

ASlope 

KS(E-Sp) 

All 

(5504) 

0.950 

(0.024) 

0.317 

(0.025) 

0.000 

0-3 

(5187) 

0.942 

(0.026) 

0.619 

(0.053) 

0.000 

0-0.5 

(3813) 

0.878 

(0.054) 

0.282 

(0.041) 

0.000 

0-1.45 

(2567) 

0.869 

(0.055) 

0.478 

(0.070) 

0.000 

0.5-1 

(1321) 

0.215 

(0.195) 

0.101 

(0.083) 

0.079 

1.45-3 

(2620) 

0.878 

(0.062) 

0.602 

(0.105) 

0.000 

0-0.33 

(2538) 

0.910 

(0.047) 

0.342 

(0.048) 

0.000 

0-1.2 

(1331) 

0.762 

(0.094) 

0.375 

(0.080) 

0.001 

0.33-0.66 

(2026) 

0.132 

(0.195) 

0.018 

(0.042) 

0.280 

1.2-1.8 
(2828) 

0.953 

(0.025) 

0.578 

(0.059) 

0.000 

0.66-1 

(570) 

0.145 

(0.228) 

0.113 

(0.134) 

0.120 

1.8-3 

(1028) 

0.709 

(0.171) 

0.522 

(0.181) 

0.000 

0.33-1 

(2596) 

0.332 

(0.184) 

0.075 

(0.041) 

0.089 






Figures [8] and [9] are similar to Figs. [6] and [7] respectively, 
but refer to the T-R relation for the global sample of WINGS 
galaxies (Fig. [8j) and for galaxies in two different ranges of 
the local density (Figs. [9j. Columns 6, 7 and 8 of Table [T| 
report the values of the coefficients in these cases, as well as 
the values obtained splitting the LogEio interval (0-3) into 
three parts. 

This table and the two above mentioned figures show 
that, at variance with the T -E relation, where the morpho¬ 
logical correlations turn out to be significant just in the in¬ 
ner regions of clusters, the T-R relation turns out to be 
always significant in the cluster environment (at least out to 
R ~ -R 200 ), irrespectively of the local density regime. This is 
highly suggestive that the T-R relation is more robust than 
the T -E one, and that the clustercentric distance is actually 
the driving parameter of the morphological fractions in clus¬ 
ters. The above conclusions seem to overturn the widespread 


paradigm of the Morphology-Densi ty relati on, making the 
old, a l terna tive point of view by Whitmore and Gil morel 
dl99ll. 1 19931. see also IWhitmord Il995l) to come up again. 
In Section[7]we briefly discuss this point, trying to interpret 
the strong outwards weakening of T-E as due to some mech¬ 
anism of morphological broadening/reshuffling, which could 
operate in the intermediate/outer regions of nearby clusters. 
In the meantime, it is worth mentioning that the robustness 
of our results seems to be supported by the tests and sim¬ 
ulations presented in Appendix [A] We are aware that they 
cannot fully remedy the lack of data at large clustercentric 
distances. Still, waiting for the conclusion of our program for 
very wide-field (l°xl°) imaging and spectroscopy of a rep¬ 
resentative sub-sample of the WINGS clustery, we believe 


1 We are presently gathering V-, B- and u’-band Omega- 
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Figure 9. T-R for WINGS galaxies with LogSio <1.45 (median value of the Eio distribution; left panel) and 1 .45<T.ogE i o <3 (right 
panel). The meaning of the symbols is as in Fig. [6] 


the results we present here should be a good approximation 
of those we will obtain with a more complete galaxy sample. 


6.3 T-E and T-R for different global cluster 
properties 

We now seek for possible dependencies of the T-E and T-R 
relations on global cluster properties. In particular, we in¬ 
vestigate the dependencies on the X-ray luminosity (LogTx) 
and on the velocity dispersion of galaxies in the clusters 
(Loger), both quantities being somehow linked to the total 
cluster mass. Moreover, we analyse how the presence of sub¬ 
structures in the clusters influences the strenght of the T -E 
and T-R relations. 


6.3.1 Cluster mass 


Cluster velocity dispersions were computed combining 
WINGS and literature redshifts. For all but one cluster, 
they are based on more than 20 spectroscopic members, 
with an average of 92 spectroscopic members per cluster 
llCava et al. 2009|b The X-ray luminosities (0.1 2.4 keV) 
from lEbeling et akl ll 19961 . I l 9981 . l2000l ) have been converted 
to the cosmology used in this paper. 

The WINGS clusters cover a wide range of a , typi¬ 
cally between 500 and 1100 km s _1 , and Lx, typically 0.2- 
5xl0 44 erg s _1 . 

Tables [2] and [3] are similar to Table [T] (first three rows), 
but they split the WINGS cluster sample in two different in¬ 
tervals of Log a and Log Lx , respectively. The splitting values 
are choosen trying to balance (as much as possible) both the 
width and the number of galaxies of the intervals themselves. 

In spite of the expected (sometimes large) uncertainties 
of the coefficients ACC and ASlope, Tables [2] and [3] sug¬ 
gest that the results we found for the global WINGS cluster 
sample do not depend either on the velocity dispersion, or 
on the X-ray luminosity of clusters. In fact, in both cases, 


CamOVST imaging and AAOmega+2dF@AAT spectroscopy for 
54 clusters of the WINGS original sample. 


the galaxy sub-samples obtained splitting the global WINGS 
cluster sample, again show the tendency of the T-E relation 
to hold just in the inner cluster regions, while the strength 
of the T-R relation does not seem to depend on the local 
density regime. 

To include figures similar to previous ones for all 24 
cases reported in Tables [2] and [3] would be a worthless space 
waste. Thus we decided to show, in Figs liOl and llll just some 
examples illustrating the validity of the above conclusions 
(see the figure captions for more details). 


6.3.2 Subclustering 

The T-E relation has been found by D80 to hold for both 
regular (“relaxed”) and irregular (“non relaxed”) clusters of 
the local Universe and this finding has given support to the 
idea that it is a sort of general rule in the realm of galaxies. 
Indeed, since, by definition, the T-E relation concerns the lo¬ 
cal environment, assuming such a relation to have a general 
validity implies it should be found everywhere, no matter 
which is the global structure of the cluster. On the contrary, 
the lack of T-E relation in irregular (clumpy) clusters would 
suggest that either it is just a by-product of the T-R relation 
(assuming it ever holds in such kind of clusters), or that in 
the cluster substructures (infalling groups?) the morphologi¬ 
cal segregation and processing still is on-going or to come, as 
suggested by D97 for intermediate redshift clusters (a rather 
unlikely possibility for nearby clusters, indeed). Therefore, 
it is of some interest to test, in our large sample of nearby 
clusters, the effects of subclustering on the T-E and T-R re¬ 
lations. To this aim, following the classificati on given b y the 
catal og of substructures in WINGS clusters ^RamelIa et al.l 
120071 . see Table Al therein), we have divided our cluster sam¬ 
ple in three sub-samples: (i) the sample M includes the 15 
WINGS clusters for which just the main struc ture (M) wa s 
detected by the DEDICA algorithm dPisani 1993 . 1 19961) : 
(ii) the sample S1+ includes the 40 clusters for which DED¬ 
ICA found at least one substructure (S) at the same redshift 
of the main structure; (in) the sample nd (not detected) 
contains the 22 WINGS clusters whose irregular and very 
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Figure 10. T-R relation for clusters with Logcr <2.85 (upper panels) and Logcr ^2.85 (bottom panels) in different ranges of local density: 
LogEio <1.45 (left panels) and 1.45^LogEio <3 (right panels). Even at a first glance, in both intervals of Logcr the T-R relation appears 
to hold over the whole range of local density. The meaning of the symbols is as in Fig. [61 


Table 2. Similar to Table [T] but for two different intervals of Logcr 


T-E T-R 


Logcr 

Rb 

ACC 

ASlope 

KS(E-Sp) 

LogEio 

ACC 

ASlope 

KS(E-Sp) 


All 

(2129) 

0.927 

(0.034) 

0.368 

(0.040) 

0.000 

0-3 

(2112) 

0.937 

(0.029) 

0.593 

(0.054) 

0.000 

<2.85 

0-0.5 

(1222) 

0.757 

(0.113) 

0.277 

(0.064) 

0.000 

0-1.45 

(1334) 

0.871 

(0.054) 

0.527 

(0.074) 

0.000 


0.5-1 

(816) 

0.289 

(0.226) 

0.100 

(0.077) 

0.425 

1.45-3 

(778) 

0.854 

(0.096) 

0.620 

(0.153) 

0.000 


All 

(3106) 

0.862 

(0.066) 

0.298 

(0.042) 

0.000 

0-3 

(3075) 

0.936 

(0.031) 

0.690 

(0.077) 

0.000 

^2.85 

0-0.5 

(2591) 

0.812 

(0.085) 

0.280 

(0.054) 

0.000 

0-1.45 

(1233) 

0.746 

(0.113) 

0.430 

(0.117) 

0.000 


0.5-1 

(505) 

0.083 

(0.192) 

0.075 

(0.122) 

0.083 

1.45-3 

(1842) 

0.961 

(0.022) 

0.727 

(0.073) 

0.000 
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Figure 11. T-S relation for clusters with Log Lx <44.15 (upper panels) and LogTx >44.15 (bottom panels) in different ranges of 
cluster-centric distance: Rg <0.5 (left panels) and 0.5< Rg <1 (right panels). Again, the T-S relation appears to hold just in the inner 
cluster regions, irrespective of the cluster X-ray luminosity. The meaning of the symbols is as in Fig. [6] 


Table 3. Similar to TablejT) but for two different intervals of Log (Lx) 


T-S T-R 


Log(Lx) 

Rg 

ACC 

ASlope 

KS(E-Sp) 

LogS io 

ACC 

ASlope 

KS(E-Sp) 


All 

(2363) 

0.883 

(0.054) 

0.315 

(0.039) 

0.000 

0-3 

(2192) 

0.887 

(0.049) 

0.635 

(0.077) 

0.000 

<44.15 

0-0.5 

(1524) 

0.710 

(0.117) 

0.256 

(0.068) 

0.000 

0-1.45 

(1223) 

0.852 

(0.062) 

0.550 

(0.090) 

0.000 


0.5-1 

(628) 

0.276 

(0.223) 

0.136 

(0.108) 

0.236 

1.45-3 

(969) 

0.885 

(0.075) 

0.949 

(0.176) 

0.000 


All 

(3041) 

0.943 

(0.025) 

0.297 

(0.027) 

0.000 

0-3 

(2895) 

0.936 

(0.031) 

0.615 

(0.067) 

0.000 

>44.15 

0-0.5 

(2228) 

0.878 

(0.054) 

0.258 

(0.038) 

0.000 

0-1.45 

(1277) 

0.740 

(0.116) 

0.432 

(0.107) 

0.000 


0.5-1 

(653) 

0.060 

(0.225) 

0.042 

(0.120) 

0.122 

1.45-3 

(1618) 

0.944 

(0.033) 

0.562 

(0.070) 

0.000 
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M S1+ nd M S1+ nd 
Cluster samples 


Figure 12. Coefficients ASlope (upper panels) and ACC (lower 
panels) of the T-E (left panels) and T-R (right panels) relations 
for the three cluster subsamples ( M, S1+ and nd) defined in 
Section 16.3.21 The full dots and squares (red and blue, in the 
electronic version) refer to different ranges of Rg and Log Eio, for 
the T -E and T-R relations, respectively. The dashed and dotted 
lines (again red and blue, in the electronic version) correspond 
to the values found in the different cases for the whole sample of 
WINGS clusters. 

clumpy structure prevented DEDICA from detecting any 
significant structure inside them. For these cluster samples, 
Figure [12] reports the coefficients ASlope (upper panels) and 
ACC (lower panels) of the T-E (left panels) and T-R (right 
panels) relations, in two ranges of Rb and Log Eio, respec¬ 
tively. The dashed and dotted lines (red and blue, in the 
electronic version) correspond to the values found for the 
whole sample of WINGS clusters in the different cases. 

Some interesting things come to light from Figure 1121 
(a) in the inner part of clusters (Rb <0.33) and in the re¬ 
gions of high local density (LogEio >1.45), the T-E and 
T-R relations (respectively) turn out to be quite strong in 
all three sub-samples (full dots in the figure; red in the elec¬ 
tronic version). In particular, in each case the strength of the 
relations is similar to that of the corresponding one relative 
to the whole cluster sample, no matter which is the subclus¬ 
tering level (dashed lines; see also Tablc[T]|; (b) for Rg >0.33, 
the relation T-E turns out to be as strong as in the inner 
part of clusters just for the subsample M, which presumably 
contains the most regular (“relaxed”) clusters, while for the 
other two sub-samples (S1+ and nd) the lack of relation 
is confirmed at a significance level even greater than in the 
case of the whole cluster sample; (c) in the lower bin of local 
density (LogEio <1.45), the T-R relation remains strong for 
both the M and S1+ sub-samples, becoming weaker just for 
the sub-sample nd. 

In Figure [TiH the T-E and T-R relations (upper and bot¬ 
tom panels, respectively) for the regular (left panels; sample 
M) and very irregular (right panels; sample nd) clusters are 
illustrated for some particularly interesting cases. The up¬ 
per panels clearly show that: (i) for regular clusters the T-E 
relation holds outside the very inner regions too (at variance 
with our finding relative to the global cluster sample); (ii) 
in the inner cluster regions the T-E relation holds even for 
very irregular clusters (not a trivial thing, indeed). The bot¬ 


tom panels show that the T-R relation, in the whole range 
of Eio, holds for both regular and very irregular clusters. 

6.4 T-E and T-R for different stellar masses 

IPeng et all (120101). using the surveys S PSS and zCOSMOS 

1 Scoville et al.l 120071 : iLillv et all I2007I ). have analysed the 

fraction of quenched galaxies in the general field as a func¬ 
tion of both the stellar mass and the environment, over the 
cosmic time. More closely in connection with our investiga¬ 
tion, the dependence of global morphological fractions on 
galaxy stell ar mass in clusters at different redshifts has been 
analysed bv lVulcani et al.i (1201 ll ). However, the possible de¬ 
pendence of the T-E and T-R relations on M* has not been 
investigated so far, a part from some particu lar aspects of the 
problem analysed bv lBamford et all (I2009I ). In this analysis 
the fraction of early-type galaxies is found to increase with 
the local density in each bin of stellar mass, the slope of the 
Fraction-LogEio relation being almost the same in all bins. 
However, this study is based on the morphol ogical classifi¬ 
catio ns provided by the Galaxy Zoo project dLintott et al.l 
B), being thus limited by the lack of distiction between 
ellipticals and SO galaxies._ 

Stellar masses of WINGS galaxies dFritz et al.l l201ll) 
have been determined by fitting the optical spectrum (in the 
range ~3600-^70 00A) with the spe ctro-photometric model 
fully described in iFritz et aid d2007ll . In this model, all the 
main spectro-photometric features are reproduced by sum¬ 
ming the theoretical spectra of Simple Stellar Population 
(SSP) of 13 different ages (from 3xl0 6 to ~14xl0 9 years) 
and assuming a Salpeter IMF. Dust extinction is allowed 
to vary as a function of SSP age and the metallicity can 
vary among three values: Z=0.004, Z=0.02 and Z=0.05. 
These mass estimates were then corrected for color gra¬ 
dients within each galaxy (i.e. for the difference in color 
within the fibr e and over a l Okpc d iameter) and were con¬ 
verted to the iKroupal d200ll ) IMF. IFritz et aid d201ll) pro¬ 
vide mass estimates for ~5,300 WINGS galaxies, 1,540 of 
them belonging t o the galaxy sample we use in this paper. 
In lYulcani et al.l d201ll ) it is shown that these mass esti¬ 
mates are in fairly good agreement with those derived us- 
ing the rest-frame (B — V) colors and the recipe given in 
iBell and de Jonal (1200 ll ). 

Tablc[4]is similar to previous Tables[2]and[3] but it splits 
the WINGS galaxies in two intervals of Log(M*). In this case 
the splitting value roughly coincides with the median of the 
Log(M*) distribution. Provided that the T-E for this re¬ 
duced galaxy sample proves again itself to be weak or absent 
outside the inner cluster regions, from Table [I] it turns out 
that the strength of the T-E relation is greater in the high- 
mass than in the low-mass bin. We decided to further explore 
this point dividing the stellar mass interval in four relatively 
narrow bins of size 0.5, starting from Log(M*)=10 0 . Fig¬ 
ure [TT] illustrates the results. It clearly shows that, while for 

2 The spectroscopic magnitude limit of the WINGS survey 
is V=20, co r respon ding to a mass limit of Log(M*)=9.8. 
lYulcani et alj d20Ilf) actually report that the completeness of 
mass limited samples in WINGS is reached for Log(M*)=10.5. 
Thus, the lowest mass bin in Fig. 1 14l could still suffer from some 
bias. 
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Figure 13. T -E and T-R relations (upper and bottom panels, respectively) for regular (left panels; sample M) and very irregular (right 
panels; sample nd) clusters. The Rb and Eio intervals tried out in each case are indicated on the top of each panel. The meaning of the 
symbols is as in Fig. [6] 


Table 4. Similar to Table [T] but for two different intervals of Log(M*) 


T-E T-R 


Log (M*) 

Rb 

ACC 

ASlope 

KS(E-Sp) 

LogEio 

ACC 

ASlope 

KS(E-Sp) 


All 

(716) 

0.508 

(0.185) 

0.195 

(0.096) 

0.001 

0-3 

(716) 

0.697 

(0.149) 

0.510 

(0.149) 

0.000 

<10.76 

0-0.5 

(523) 

0.679 

(0.150) 

0.290 

(0.103) 

0.009 

0-1.45 

(365) 

0.467 

(0.213) 

0.404 

(0.240) 

0.006 


0.5-1 

(182) 

-0.020 

(0.337) 

-0.057 

(0.327) 

0.168 

1.45-3 

(351) 

0.638 

(0.197) 

0.554 

(0.210) 

0.002 


All 

(824) 

0.879 

(0.055) 

0.354 

(0.053) 

0.000 

0-3 

(809) 

0.876 

(0.071) 

0.681 

(0.118) 

0.000 

>10.76 

0-0.5 

(590) 

0.738 

(0.120) 

0.280 

(0.081) 

0.004 

0-1.45 

(395) 

0.753 

(0.134) 

0.644 

(0.207) 

0.001 


0.5-1 

(225) 

-0.183 

(0.289) 

-0.096 

(0.186) 

0.404 

1.45-3 

(414) 

0.772 

(0.142) 

0.637 

(0.190) 

0.001 
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Figure 14. ASlope (upper panel) and ACC (lower panel) as a 
function of Log(M*) for both T -S (full dots; red in the electronic 
version) and T-R (full squares; blue in the electronic version) 


the T-R relation both ASlope and ACC turn out to be al¬ 
most stable in the four bins of Log(M*), they are strongly 
increasing functions of the galaxy stellar mass for the T -E 
relation. In particular, in the lowest mass bin, both ASlope 
and ACC indicate that the T -E is almost absent. 

We have also tried to account for spectroscopic incom¬ 
pleteness by applying a statistical correction to our sample. 
This is obtained by weighting each galaxy by the inverse 
of the ratio of the number of spectra yielding a redshift to 
the total number of galaxies in th e photometric catalogue, 
in bins of 1 mag rtCava et al.l 120091) . We do not report here 
the results of such additional analysis. We just mention that 
they turn out to be quite similar to those shown in Table 0] 
and in Figure fl4l 

This might induce to draw the conclusion that the lack 
of T -E we observe outside the inner cluster regions is due 
to an excess in these regions of galaxies in the lowest stel¬ 
lar mass bin with respect to the other bins. However, this 
conclusion turns out to be ruled out by the 2S-KS applied 
to the clustercentric distances of galaxies in the four mass 
bins. _ 

It is well known (see for instance IVulcani et al.l 1201 ll ) 
that the three broad morphological types (E/SO/Sp) have 
quite different stellar mass distributions. This is confirmed 
by the 2S-KS applied to our galaxy sample. It might be spec¬ 
ulated that this fact, combined with t he dependence of th e 
mass distribution on the local density (IVulcani et al.ll2012l) . 
is fully responsible of the very existence of the T -E relation. 
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Figure 15. T-E (upper panel) and T-R (lower panel) for galaxies 
in the Log(M*) bins (11—11.3) and (10.3-10.6), respectively. 


In other words, the T -E could just be a consequence of the 
combined dependence of both the morphology and the local 
density on the stellar mass. This hypothesis would obviously 
imply that the T-E relation should not be observed at any 
(fixed) galaxy mass. 

The upper panel of Figure [15] contradicts such an ex¬ 
pectation. It shows that, in spite of the poor statistics (218 
galaxies), the T -E turns out to be remarkably strong even 
in quite narrow stellar mass bins [ll<Log(M*)<11.3]. For 
the sake of clarity, in this plot we use the histograms to 
represent the T-E, instead of the points, since the big error 
bars associated to the points whould weight down too much 
the figure. It is also worth noticing that we have chosen the 
narrow mass bin in the high-mass part of the mass distribu¬ 
tion, since the T-E progressively weakens towards low-mass 
galaxies (see Figure [Til) . The lower panel of Figure [15] shows 
that, similarly to T-E, the T-R relation appears to be strong 
enough also in a narrow stellar mass bin, while, at variance 
with the T-E (and according to Figure [T4]l . it is quite strong 
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also for galaxies in the low-mass region of the distribution 
[10.3<Log(M*)<10.6]. 


7 SUMMARY AND OPEN ISSUES 

In the present analysis, performed using the WINGS 
database of galaxies in nearby clusters, we have shown that: 

• the correlation between morphological fractions and lo¬ 
cal density (T-E) exists only in the very inner regions of 
nearby clusters, almost vanishing outside I/ 3 XR 200 , apart 
from very regular (non substructured) clusters, for which 
the T-E relation also holds outside the inner regions; 

• in contrast, the strenght of the correlation between mor¬ 
phological fractions and clustercentric distance ( T-R ) re¬ 
mains almost unchanged over the whole range of local den¬ 
sity for both regular and irregular clusters, only slightly less¬ 
ening for extremely clumpy clusters (bad center determina¬ 
tion?); 

• a couple of suitable tests and two different sets of nu¬ 
merical simulations support the validity of these results even 
considering the possible biases arising from the limited clus¬ 
ter area coverage of the WINGS imaging; 

• the above findings hold irrespective of both the global 
cluster mass (velocity dispersion and X-ray emission) and 
the stellar mass of galaxies; 

• the strength of the T-E relation (where present) in¬ 
creases with increasing galaxy stellar mass, while this effect 
is not found for the T-R relation; 

• both the T-E (where present) and the T-R relations are 
remarkably strong even in quite narrow stellar mass bins. In 
particular, for the T-E relation, this rules out the hypoth¬ 
esis that it could just be a consequence of the combined 
dependence of both morphology and local density on stellar 
mass. 


These results could lead us to conclude that the param¬ 
eter actually driving morphological fractions in nearby clus¬ 
ters is the distance from the cluster center, rather than the 
local density, as commonly believed. In this scenario the T-E 
relation would just be a by -product of the T-R relation, as 
already claimed bv I Whitmore and Gilmore ( 1993 1. Against 
this conclusion, how ever , one can br ing forward the argu¬ 
ment (challenged bv I Whitmore|[l995h that the T-E relation 
is also found in the general field (groups+pairs+single galax¬ 


ies; Bhavsaill98ll;[de Souza et al. ll 19821 ; IPostman and Geilei] 


I 1984k Helsdon and Ponman 20031 '!. In fact, according to the 
standard paradigm, galaxy clusters are progressively built 
up through continuous infalling into the main cluster halo 
of galaxies coming from the surrounding field. Thus, it is 
worth asking oneself whether the lack of the T-E in the out- 
of-center part of clusters really implies the fall of the T-E 
paradigm or, instead, some other mechanism should be in¬ 
voked to interpret our findings. 

First, we note that the results illustrated in this paper 
(Section 16.3.21 in particular) suggest that the T-E relation 
holds only in dynamically evolved regions of nearby clusters, 
i.e. the whole clusters, or just their inner parts, for regular 
(“relaxed”) or substructured (“non relaxed”) clusters, re¬ 
spectively. Therefore, it seems to us that the apt question 
is: what is causing the observed T-E weakening in dynami¬ 
cally “non relaxed” regions? We speculate that some mecha¬ 


nism of morphological broadening/redistribution is respon¬ 
sible for the T-E weakening in the intermediate regions of 
substructured (“non relaxed”) clusters. 

Two different (and perhaps complementary) ways to at¬ 
tain this effect are conceivable: (i) general field galaxies lose 
(release) the dependence of their morphology on the local 
density when they infall into the cluster, recovering it only 
when (and where) the dynamical equilibrium has been al¬ 
ready reached (only inner regions for substructured, “non 
relaxed” clusters); (ii) galaxies having already experienced 
the very dense interiors of clusters and moving back to the 
intermediate/outer regions, or ’free’ galaxies coming in these 
regions from different sub-halos are ‘T-E untied’, thus dilut¬ 
ing the relation likely existing in the infalling galaxy popu¬ 
lation. 

The first mechanism could result from reshuffling of 
galaxy locations due the gravitational interactions within 
the cluster. The second mechanism envisages a scenario 
in which galaxies responsible for morphological broadening 
could come from the innermost cluster regions (likely af¬ 
ter one or more passages close to the cluster center) and/or 
could be detached from their original infalling halos, there¬ 
fore losing their pristine correlation between morphology 
and E. 

Of course, the two outlined mechanisms of morpho¬ 
logical broadening could operate in tandem. Unfortunately, 
there are no observational hints about which one of them is 
prevailing, or even about their very existence. As a matter 
of fact, it is presently unknown which process is responsible 
for the weakening of T-S in the intermediate cluster regions. 
The increasingly powerful and sophisticated A CDM mod¬ 
els of galaxy formation and evolution inside DM halos could 
help shed some light on these findings. 
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APPENDIX A: TESTING THE EFFECTS OF 
THE LIMITED CLUSTER AREA COVERAGE 

The following three tests heve been devised with the aim of 
investigating the effects of the irregular and/or limited clus¬ 
ter coverage of the WINGS imaging on the results we present 
about the T -S relation in Section [6T] and, marginally, about 
the T-R in Sec. 16.21 

Al Irregular Image Shape 

In the first test, using the same galaxy sample of Sec. [6l 
we produce a new versions of the T -E relation in which 
each galaxy, with its proper clustercentric distance (Rb), 
is weighted according to the inverse of its circumferential 


Table Al. Similar to Table [T] but for just T-T and for two 
different values of the minimum coverage fraction (CF) 


T-T. 

CF 

Rb 

ACC 

ASlope 

KS(E-Sp) 


0-1 

(1710) 

0.853 

(0.065) 

0.352 

(0.057) 

0.000 

CF=0.5 

0-0.5 

(3351) 

0.878 

(0.056) 

0.282 

(0.041) 

0.000 


0.5-1 

(454) 

-0.087 

(0.278) 

0.022 

(0.123) 

0.210 


0-1 

(381) 

0.782 

(0.106) 

0.512 

(0.137) 

0.000 

CF=0.75 

0-0.5 

(2551) 

0.885 

(0.051) 

0.286 

(0.040) 

0.000 


0.5-1 

(127) 

-0.327 

(0.290) 

-0.138 

(0.394) 

0.037 


coverage, that is the fraction of circumference (of radius Rb) 
covered by the image. This is equivalent to assume that, for 
each galaxy falling into the image, other (similar) galaxies 
could exist at the same clustercentric distance, but falling 
outside the image. This weighting also (obviously) increases 
the ’nominal’ number of galaxies (in particular spirals), part 
of them being actually virtual objects. This test tries to 
partially account for the non circular (and irregular) shape 
of our images (from INT, in particular). 

Figure 1X11 is similar to Fig. [7] but it is obtained using 
the above outlined weighting procedure of galaxies. From 
the comparison between figures IT] and [All we conclude that 
the non circular and irregular shape of the images should 
not invalidate the conclusion reported in Sec. that the 
strenght of the T-T, relation in clusters strongly depends on 
the clustercentric distance. 

A2 Limited cluster area coverage 

In the second test, we compare the T-T obtained using 
galaxies in the whole cluster sample with the corresponding 
ones in which only galaxies in clusters with large values of 
the cluster area coverage are included. This is equivalent to 
bias the analysis towards smaller and more distant clusters. 

Table [Al] is similar to Table |TJ but it reports the re¬ 
sults for just the T-T and for two different values of the 
minimum allowed coverage fraction (CF>0.5 and CF>0.75). 
More precisely, to assume for instance CF>0.5, means that 
just galaxies in clusters whose CF in the considered range of 
Rb is greater than 0.5 are each time included in the sample. 
This obviously reduces the number of employed galaxies, but 
makes the T-T relation more robust against the CF issue. 
From the comparison between Tables [II and I A1 1 we conclude 
again that the limited cluster area coverage of the images 
does not invalidate the conclusion reported in Sec. 16.21 that 
the strenght of the T -E relation in clusters strongly depends 
on the clustercentric distance. 

A3 Numerical Simulations 

Finally, we use two different sets of numerical simulations to 
investigate the effect of the limited cluster area coverage on 
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Log (£,„) Log (S 10 ) 


Figure Al. T-S for WINGS galaxies with Rg <0.5 (in units of R 200 ; left panel) and 0.5 Rg <1 (right panel). The galaxies have been 
weighted according to the inverse of their circumferential coverage (see text). The meaning of the symbols is as in Fig. [6] 


T-£ (0.5SR b <1) 



ACC ASlope 


Figure A2. Results of the simulations SSI (upper panels) and 
SS2 (lower panels) for the T- S in the R interval (0.5-1). The left 
and right panels illustrate the distributions of ACC and ASlope, 
respectively. The dark and light histograms refer to the HR and 
HEjo hypotheses, respectively. 


the T- E in different ranges of the clustercentric distance. In 
order to perforin the simulations we have to assume some 
dependence of the morphological fractions on either R or 
E 10 (or even on both). Therefore, these simulations can also 
be used, in principle, to test which one of the two (or three) 
hypotheses is able to reproduce the observed T-E and T-R 
relations. 

The first set (Simulation Set 1: SSI) totally ignores the 


physics of the clusters. It simply produces artificial samples 
of points with random polar angles (circular symmetry) and 
radial coordinates randomly extracted from a distribution 
suitably chosen in order to reproduce the Rg-LogEio rela¬ 
tion observed in the WINGS galaxy sample. Once the E 10 
values of the simulated points have been computed in the 
usual way, the artificial sample is randomly thined out fol¬ 
lowing a depletion law depending on the radial coordinate 
and suitably chosen in order to reproduce the variation of 
the total number density of the WINGS sample as a func¬ 
tion of Rg. This procedure should properly account for the 
incompleteness of the real sample due to both the irregular 
image shape and the limited coverage fraction. However, the 
lack of physics, in particular of gravitational clustering ( 2 - 
point or 3-point correlation function), makes the scatter of 
E 10 in the simulated R-LogEio relations lower than in the 
real sample of WINGS galaxies. We compensate this lack of 
physics through an additional scatter on E 10 , tuned on the 
observed R LogEio relation and depending on the radial co¬ 
ordinate. Then, we consider the two alternative hypothesis 
that the morphological fractions (Fe and Fs) depend either 
just on E 10 (Hypothesis E 10 : HE 10 ) or just on R (HR). We 
assume in the two cases the probability that the points in 
our simulated samples were elliptical or spiral galaxies to co¬ 
incide with the weighted linear fits of the relations in Fig. [ 6 ] 
[.Fe(Eio) and Rs(Eio)] or in Fig. [ 8 ] [F e (Rb) and Fs(Rg)], 
respectively. They are: 

F e = 0.125 + 0.136 L 03 E 10 ; F s = 0.493 - 0.177 LogEio 

F e = 0.401 - 0.25 Rg; F s = 0.091 + 0.368 Rg 

We also formally contemplate the ’mixed’ hypothesis 
(HMIX), in which the morphological fractions depend on 
both E 10 and R. In this case we roughly assume: 

F e = [F e { E 10 ) + F E (R)}/ 2; F s = [F S (T, W ) + F S (R)}/ 2 . 

In the second set of simulations (SS2) we adopt the 
same procedure to randomly assign the morphological type, 
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but we use the real sample of WINGS galaxies to take the 
(measured) values of Eio and R. 

For both kinds of simulations (SSI and SS2) and for 
each hypothesis about the dependence of the morphological 
fractions (HEio, Hi?., HMIX), we make several throws, each 
throw producing T -E and T-R relations similar to those il¬ 
lustrated in Figures I6I7I8I9I (R intervals: 0-1, 0-0.5, 0.5-1; 
LogEio intervals: 0-3, 0-1.45, 1.45-3), with the same num¬ 
ber of points (galaxies). Moreover, for each throw we record 
the values of the coefficients ACC and ASlope defined in 
Sec. 16.11 Finally, in each case, we compare the distributions 
of ACC and ASlope with the corresponding values obtained 
for the real sample of WINGS galaxies (see the previously 
mentioned figures). In this way, for both SSI and SS2, we 
can associate to each hypothesis (HEio, Hi?., HMIX) and 
for both the T -E and T-R relations (in the corresponding 
intervals of R and Eio, respectively) a set of deviates of the 
observed coefficients ACC and ASlope (CCdev and SLdev; 
both in r.m.s. units) with respect to the corresponding dis¬ 
tributions, as well as a set of probabilities (Pec and P sl) 
that these coefficients are extracted from the distributions 
themselves. 

Tables IA2I and IA3I illustrate the results obtained for 
the two sets of simulations (SSI and SS2, respectively). For 
each hypothesis (HEio, Hi?, HMIX), the tables report the 
deviates and the probabilities of the coefficients ACC and 
ASlope of the observed T -E and T-R relations for differ¬ 
ent intervals of i? and LogEio, respectively. They are com¬ 
puted comparing each coefficient with the corresponding, 
proper distribution obtained from the simulations. It should 
be noted that, due to the rather cumbersome (and CPU time 
consuming) procedure needed to produce a single SSI sam¬ 
ple , the number of throws used to obtain the distributions 
of coefficients CCdev, Pec, SLdev and Psi, is almost one 
order of magnitude smaller for SSI (N s i m =100-i-200) than 
for SS2 (N s i m =1000). This fact, together with the use of the 
real WINGS quantities (Eio and Rb ) in SS2, suggests that 
the robustness of the results is greater for SS2 than for SSI. 
In both Table [A~2l and Table [Ail we report in boldface those 
values of CCdev and SLdev larger than 3, as well as those 
values of Pec and P sl lower that 0.05. Moreover, we mark 
with one or two asterisks those values of Pec and P sl equal 
to or less than 0.01, respectively. 

In spite of the quite rough procedure used to generate 
the simulated samples in the SSI and to attribute the mor¬ 
phological type to each simulated point (in both SSI and 
SS2), Tables IA2I and IA3I provide us with some indication 
about the driving parameter of the morphological fractions 
in clusters: 

(%) the ’canonical’ hyphotesis that the morphological 
fractions in clusters just depend on the local density (HEio) 
is not able to reproduce the observed T-R relation, either 
for the whole sample, and for the two Eio intervals. This is 
particularly evident in Table [~\3l 

(ii) although specifically built to obey the observed 
morphology-density relation, the HEio simulations just 
marginally reproduce the very T -E relation in the outer 
part of clusters (i? interval: 0.5-1), usually producing much 
steeper and more correlated distributions than the observed 
one (see light histograms in Figure |X2j. While for SSI this 
might be thougth as a consequence of the quite regular shape 
of the simulated clusters (see Section 16.3.21) . such an ex¬ 


planation is ruled out by SS2, where the observed WINGS 
quantities (Eio and Rb) have been used; 

(in) on the contrary, the alternative hyphotesis that the 
morphological fractions in clusters just depend on the clus¬ 
tercentric distance (HJ?) quite well reproduces both the T-R 
relation (in any Eio interval) and the T-E relations in the 
outer cluster regions (see dark histograms in Figure IX2l) . but 
it seems to fail in reproducing the T-E in the inner cluster 
regions. This is particularly true for the SSI (Table lA2l) ; 

(iv) the mixed hypothesis that the morphological frac¬ 
tions actually depend on both Eio and R (HMIX) turns out 
to be fairly consistent with all observed T-E and T-R re¬ 
lations, but (perhaps) the T-R in the full LogEio interval 
(0-3). 

The above remarks seem to favour the HJ? hypothesis 
with respect to the classical HEio scenario, although this re¬ 
sult might critically depend on the simulation procedure and 
might be modified, for instance, by a different (not linear) 
representation of the relations Tb(Eio), Ts(Eio), Fe{R ) and 

F S (R). 
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Table A2. Simulation Set 1 (SSI). Since the simulations of this set turn out to be rather consuming in terms of CPU time, the number 
of throws used here to obtain the distributions of coefficients CCdev, P cCi SLdev and Ps£ varies from 100 to 200. 


T- S T-R 


HEio 


Rb interval 

CCdev 

P CC 

SLdev 

P SL 

LogEio interval 

CCdev 

P CC 

SLdev 

P SL 

0-1 

1.13 

0.32 

0.10 

0.90 

0-3 

1.10 

0.14 

4.02 

0.00** 

0 - 0.5 

1.50 

0.24 

0.61 

0.49 

0 - 1.45 

1.42 

0.08 

1.90 

0.07 

0.5 - 1 

4.04 

0.01* 

1.85 

0.03 

1.45 - 3 

1.11 

0.12 

1.82 

0.02 

HR 

Rb interval 

CCdev 

p CC 

SLdev 

P SL 

LogEio interval 

CCdev 

P CC 

SLdev 

P SL 

0-1 

0.11 

0.86 

1.79 

0.02 

0-3 

1.18 

0.22 

0.04 

0.92 

0 - 0.5 

0.83 

0.23 

3.21 

0.00** 

0 - 1.45 

0.12 

0.92 

0.87 

0.34 

0.5 - 1 

0.00 

0.99 

0.50 

0.65 

1.45 - 3 

0.17 

0.82 

0.07 

0.94 

HMIX 

Rb interval 

CCdev 

P CC 

SLdev 

Psx 

LogEio interval 

CCdev 

P CC 

SLdev 

Psx 

0-1 

0.10 

0.96 

0.93 

0.26 

0-3 

0.29 

0.66 

1.56 

0.10 

0 - 0.5 

0.14 

0.93 

0.99 

0.33 

0 - 1.45 

0.70 

0.32 

0.47 

0.57 

0.5 - 1 

1.44 

0.31 

0.61 

0.60 

1.45 - 3 

0.69 

0.36 

1.29 

0.25 


Table A3. Simulation Set 2 (SS2). The simulations of this set are much faster than in the case of SSI. Therefore, the number of throws 
used here to obtain the distributions of coefficients CCdev, Pecs SLdev and P sl is everywhere 1000. 


T -E T-R 


HS 10 


Rb interval 

CCdev 

P CC 

SLdev 

Psl 

LogEio interval 

CCdev 

P CC 

SLdev 

P SL 

0-1 

1.11 

0.33 

0.19 

0.86 

0-3 

1.14 

0.18 

4.31 

0.00** 

0 - 0.5 

1.47 

0.24 

0.51 

0.61 

0 - 1.45 

1.33 

0.07 

3.24 

0.00** 

0.5 - 1 

5.64 

0.00** 

2.16 

0.04 

1.45 - 3 

1.21 

0.11 

3.23 

0.00** 

HR 

Rb interval 

CCdev 

p CC 

SLdev 

P SL 

LogEio interval 

CCdev 

P CC 

SLdev 

P SL 

0-1 

0.28 

0.76 

2.29 

0.02 

0-3 

1.72 

0.12 

0.05 

0.94 

0 - 0.5 

0.66 

0.44 

2.72 

0.01* 

0 - 1.45 

1.72 

0.18 

1.53 

0.14 

0.5 - 1 

0.03 

0.97 

0.51 

0.61 

1.45 - 3 

0.98 

0.41 

0.13 

0.92 

HMIX 

Rb interval 

CCdev 

P CC 

SLdev 

P.SL 

LogEio interval 

CCdev 

P CC 

SLdev 

P SL 

0-1 

0.30 

0.77 

1.25 

0.18 

0-3 

0.00 

0.99 

2.21 

0.03 

0 - 0.5 

0.21 

0.82 

1.19 

0.26 

0 - 1.45 

0.34 

0.73 

0.84 

0.41 

0.5 - 1 

1.76 

0.18 

0.74 

0.44 

1.45 - 3 

0.48 

0.61 

1.59 

0.10 


© 2012 RAS, MNRAS 000, [T] ?? 




























